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Abstract 

We propose an alternative to the usual time-independent Born-Oppenheimer ap- 
proximation that is specifically designed to describe molecules with symmetrical Hy- 
drogen bonds. In our approach, the masses of the Hydrogen nuclei are scaled differently 
from those of the heavier nuclei, and we employ a specialized form for the electron en- 
ergy level surface. Consequently, anharmonic effects play a role in the leading order 
calculations of vibrational levels. 

Although we develop a general theory, our analysis is motivated by an examination 
of symmetric bihalide ions, such as FHF~ or ClHCl~ . We describe our approach for 
the FHF~ ion in detail. 
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1 Introduction 



In standard Born-Oppenheimer approximations, the masses of the electrons are held fixed, 
and the masses of the nuclei are all assumed to be proportional to e~ 4 . Approximate solutions 
to the molecular Schrodinger equation are then sought as expansions in powers of e. For 
the time-independent problem, the electron energy level surface is also assumed to behave 
asymptotically like a quadratic function of the nuclear variables near a local minimum. 

In this paper and in a future one jlj, we propose an alternative approximation for 
molecules that contain Hydrogen atoms as well as some heavier atoms, such as Carbon, 
Nitrogen, or Oxygen. Our motivation is to develop an approach that is specifically tailored 
to describe the phenomenon of Hydrogen bonding. 

In this paper, we examine the specific case of systems with symmetric Hydrogen bonds, 
such as FHF~ . In [4 , we plan to study non-symmetric cases, where the structure of the 
typical electron energy surface is very different. The mathematical analysis of that situation 
is consequently completely different. 

The model we present here differs from the usual Born-Oppenhimer model in two ways: 

1. We scale the masses of the Hydrogen nuclei as e~ 3 instead of e~ 4 . This is physically 
appropriate. If the mass of an electron is 1, and we define e -4 to be the mass of a C 12 
nucleus, then e = 0.0821, and the mass of a H 1 nucleus is 1.015 e~ 3 . 

2. We do not assume that the electron energy level is well approximated by an e-independent 
quadratic function near a local minimum. Instead, we allow it to depend on e and to 
take a particular form that we specify below. The particular form we have chosen is 
motivated by a detailed examination of the lowest electronic potential energy surfaces 
for FHF~ and ClHCl' . 

Although symmetric bihalide ions are quite special, our approach is flexible enough to 
describe more general phenomena. For example, the lowest electron energy surface for FHF~ 
has a single minimum with the Hydrogen nucleus mid-way between the two Fluorines. Our 
model can handle situations with single or double wells in the coordinates for a Hydrogen 
nucleus that paticipates in Hydrogen bonding. We hope that the ideas in this paper and jlj 
might provide some insight into some properties of Hydrogen bonded systems. 

Our model leads to a different expansion from the usual Born-Oppenheimer approxima- 
tion. For Hydrogen nuclei not involved in Hydrogen bonding, the vibrational energies are 
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of order e 3//2 , while the vibrational energies for the other nuclei and the Hydrogen nuclei in- 
volved in the symmetric Hydrogen bonding are of order e 2 . Furthermore, anharmonic effects 
must be taken into account for a Hydrogen nucleus involved in Hydrogen bonding at their 
leading order, e 2 . In the standard Born-Oppenheimer model, all vibrational energies appear 
in a harmonic approximation at order e 2 . Anharmonic corrections enter at order e 4 . 

We present our ideas only in the simplest possible situation. In that situation, there are 
only 3 nuclei, and they are constrained to move along a fixed line. We plan to study more 
general possibilities, such as bending of the molecule, in the future. 

The paper is organized as follows: In Section EJ we present the formal expansion. In 
Section El we state our rigorous results as Theorems 13 . 61 and 13 . 71 The proofs of some technical 
results are presented in Section HJ 

Acknowledgements The authors would like to thank Thierry Gallay for several helpful 
conversations. George Hagedorn would like to thank T. Daniel Crawford for teaching him 
to use the Gaussian software of computational chemistry. 

2 Description of the Model 

We study a triatomic system with two identical heavy nuclei A and B, and one light 
(Hydrogen) nucleus C. We begin by describing the Hamiltonian for this system in Ja- 
cobi coordinates. We let xa and xb he the positions of the heavy nuclei, and let xc be 
the position of the light nucleus C. We let their masses be tjia = vn^ and mc- We 

let R = — - - — B B — — — denote the center of mass of all three nuclei, and let 

m A + m B + m c 

xab = %A ^ XB denote the center of mass of the heavy nuclei. We let W = xb — %a 
be the vector from nucleus A to nucleus B and let Z = xq — %ab be the vector from the 
center of mass of A and B to C. We assume the electronic Hamiltonian h e only depends on 
the vectors between the nuclei, and we set tuab = mA + m B and M = tra + tub + mc- In 
the original variables, the Hamiltonian has the form 

1 A 1 A 1 A , , X 

A *vt - 7^ _ 7, A x c + h e (x B - x A , x c - x A , x c - x B j- 



2m A 2m B 2m c 

In these Jacobi coordinates, it has the form 

A R - n A w - A z + h e {W, Z + W/2, Z - W/2). 



2M 2 rriA ttlb 2 tuab m c 
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Since we are interested in bound states, we discard the kinetic energy of the center of mass. 
We take the electron mass to be 1, and the masses of the heavy nuclei to be ttia — tub — e~ 4 /i, 
for some fixed /i. The mass of the light nucleus is m c = e~ 3 u, for some fixed v. The electronic 
Hamiltonian h e then becomes h e (W, Z + W/2, Z — W/2) = h(W, Z), so that the Hamiltonian 
of interest is 

e 4 e 3 / e v \ 

- — A w - 1+ A z + h(W,Z). 

\x 2v \ 2fx J 

This computation is exact and valid in any dimension. 

To simplify the exposition, we drop the term — in the factor that multiplies Az- It 

2 \x 

gives rise to uninteresting, regular perturbation corrections. Also, for simplicity, we assume 
11 — 2 and v — 1. This can always be accomplished by trivial rescalings of W and Z. 

To describe our ideas in the simplest situation, we restrict W and Z to one dimension. 
Thus, we are not allowing rotations or bending of the molecule. Furthermore, we introduce e 
dependence of the electronic Hamiltonian to model the pecularities of symmetric Hydrogen 
bonds that we describe below. 

These considerations lead us to study the Hamiltonian 

The electron Hamiltonian h(e, W, Z) is an operator in the electronic Hilbert space that de- 
pends parametrically on (e, W, Z) and includes the nuclear repulsion terms. For convenience, 
we assume that h(e, W, Z) is a real symmetric operator. 

We now describe the specific e dependence of h(e, W, Z) that we assume. Although the 
electron Hamiltonian does not depend on nuclear masses, the parameter e is dimensionless, 
and thus may play more than one role. The dependence of h on e we allow is motivated by the 
smallness of a particular Taylor series coefficient we observed in numerical computations for 
the ground state electron energy level for the real system FHF~ . We allow only the ground 
state eigenvalue to depend on e. Otherwise, our electron Hamiltonian is e-independent. With 
the physical value of e inserted in our Hamiltonian, we obtain the true physical Hamiltonian. 

From numerical computations of E(W, Z) for FHF~, we observed that the Z 2 coefficient 
in the Taylor expansion about the minimum (W , 0) of the ground state potential energy 
surface had a small numerial value, on the order of the value of e = e , where e was defined 
by setting 4 equal to the nuclear mass of the C 12 isotope of Carbon. 
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The value of eo is roughly 0.0821. We define a 2 so that the true Z 2 Taylor series term is 
a 2 e Z 2 • We then obtain h(e, W, Z) by adding (e — e ) a 2 Z 2 to the ground state eigenvalue 
E(W, Z). We make no other alterations to the electron Hamiltonian. When e = e , our 
h(e, W, Z) equals the true physical electron Hamiltonian h(e , W, Z). 

Thus, we assume the ground state electron level has the specific form 

E 1 (e, W, Z) = E + a^W-Wo) 2 + ( a 2 e - a 3 (W - W ) ) Z 2 + a 4 Z 4 + (2.2) 

with cij = 0(1). As we shall see, the leading order behavior of the energy and the wave 
functions for the molecule are determined from the terms written explicitly in ()2.2j) . The 
terms not explicitly displayed are of orders (W — Wo) a Z 2 @, where a and (3 are non-negative 
integers that satisfy a + (3 > 3. They play no role to leading order, but contribute to higher 
order corrections. 

We assume ai, a 3 , and a 4 are positive, but that a 2 can be positive, zero, or negative. 
When a 2 is negative, Ei(e, W, Z) has a closely spaced double well near (Wo, 0) instead of a 
single local minimum. 

To ensure that the leading part of Ei(e, W, Z), 

£i(e, W, Z) = E + ai {W - W ) 2 + (a 2 e - a 3 (W - W ) ) Z 2 + a 4 Z\ 
is bounded below, we assume that either 

al < 4aia 4 , (2.3) 

or 

al = AcLidi and a 2 > 0. (2.4) 

These conditions are equivalent to the property Ei(e, W, Z) > —C for some C, since we 
can write 

E^e, W,Z) = a, [{W-W Q ) - ^-Z 2 ^} + (a 4 - ^jZ A + a 2 eZ 2 . 

By rescaling with w = {W — Wo)/e and z = Z/e 1 ^ 2 , we see that the Hamiltonian 

e 4 8 2 e 3 8 2 



+ £i(e, W, Z) 



2 dW 2 2 dZ 2 

is unitarily equivalent to e 2 times the e-independent Normal Form Hamiltonian 



Id 2 Id 2 

HnF = ~ 2 d^ 2 ~ 2 M + EMW > Zl (2 ' 5) 
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where 

E NF (w, z) = ai w 2 + (a 2 - a 3 w) z 2 + a 4 z A . 



(2.6) 



Remark Although we do not use it, further scaling shows that Hnf is essentially a three- 
parameter model, since the change of variables w = as, z = at, yields 

-2 ( 1 d 2 1 2 2 2 2 4 

nf — & \ 2d? ~ 2di* +aiS + a2t ~ asSt + 1 

with 

-1/6 0,1 a 2 a 3 

a = a 4 , ax — a 2 = an d «3 " 



2/3' " z 2/3' " ^ 5/6' 

Under conditions ()2.3|) or ()2.4j) . H^p is essentially self-adjoint on C^(IR 2 ) and has purely 
discrete spectrum. This last property is easy to verify under condition (|2.3jl . or condition 
(j2.4|) with a 2 > 0, because E^p(w, z) tends to infinity as || (it?, z) || — > oo. When (j2.4|) is 
satisfied with 02 = 0, the result is more subtle because E^f(w, z) attains its minimum value 
of zero along a parabola in (w,z). In that case we prove that the spectrum is discrete in 
Proposition 13.11 

Explicit Computations for FHF~ The expression ()2.2|) is clearly special. Our compu- 
tations for FHF~ that motivate this expression have roughly the following values, where 
distances are measured in Angstroms and energies are measured in Hartrees: 

W = 2.287, 

E = -200.215, 

ai = 0.26, 

a 2 = 1.22 (if e = 0.0821), 



a 3 = 1.29, 
a 4 = 1.62. 

These results came from fitting the output from Gaussian 2003 using second order Moller- 
Plesset theory with the aug-cc-pvtz basis set. We observed that the process of fitting the 
data was numerically quite unstable, and that condition (J2.3J) was barely satisfied by these 



-0.4 
2.1 



2.2 



2.3 2.4 
F - F Distance 



2.5 



Figure 1. Contour plot of the ground state electronic potential energy surface in the Jacobi 
coordinates (W, Z). It is obviously not well approximated by a quadratic. Our technique 
exploits the flatness of the surface in the Z direction near the minimum. 

The experimentally observed values [H] for the excitation energies to the first symmet- 
ric stretching vibrational mode and the first asymmetric vibrational mode of FHF~ are 
583.05 cm -1 and 1331.15 cm" 1 , respectively. With the values of aj above, the leading order 
calculation from our model predicts 600 cm -1 and 1399 cm -1 . By leading order, we mean 
E + e 2 £2 in the expansion we present below. These values depend sensitively on precisely 
how we fit the potential energy surface, which itself depends sensitively on the electron 
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structure calculations. By comparison, Gaussian 2003 with the aug-cc-pvdz basis set pre- 
dicts harmonic frequencies of 608 cm -1 and 1117 cm -1 . We could not obtain frequencies for 
the aug-cc-pvtz basis set from Gaussian because of our computer limitations. 

For some very recent numerical results for vibrational frequencies of FHF~ that appeared 
as we were finishing this paper, see |2]. 

We now mimic the technique of [3 to obtain an expansion for the solution to the eigen- 
value problem for (|2.1|) . We could have used the technique of [S], but that would have led 
to more complicated formulas. 

For convenience, we replace the variable W by W — Wq, so that henceforth, Wq = 0. 

The technique of |3 a uses the method of multiple scales. Instead of searching directly for 
an eigenvector \l/(e, W, Z) for ()2.1j) . we first search for an eigenvector i/j(e, W, Z, w, z) for an 
operator that acts in more variables. When we have determined ip, we obtain \l/ by setting 

*(e, W, Z) = ^(e, W, Z, W/e, Z/t 1 ' 2 ). 

This is motivated physically by the following observation: The dependence of the electrons 
on the nuclear coordinates occurs on the length scale of (W, Z), while the semiclassical 
quantum fluctuations of the nuclei occur on the length scale of (w, z). To leading order in 
e, these effects behave independently. 

The equation for ip is formally 

H 2 (e) ^(e, W, Z, w, z) = £(e) V(e, W, Z, w, z), (2.7) 

where 

R(e) = _ t ^- - e 5/2 J^L _ f! ?_ 

21 ' 2 dW 2 dWdw 2 dw 2 2 dZ 2 dZ dz 2 dz 2 

+ [h(e, W, Z) - E(e, W, Z)\ + E(e, ew, e 1/2 z) 

oo 

+ £ m/2 ( T m /2{W, Z) - T m/2 (ew, e^z) ) . (2.8) 

m=6 

The functions T m / 2 in this expression will be chosen later. Different choices yield equally 
valid expansions for \l/(e, W, Z), although they alter the expressions for ip(e, W, Z, w, z) by 
converting (W, Z) dependence into (w, z) dependence. 
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In (|2.8J1 . we expand both E(e, ew, e 1 ^ 2 z) and T m / 2 (ew, e 1//2 z) in Taylor series in powers 
of e 1 / 2 . We then make the Ansatz that (|2.7j) has formal solutions of the form 

^(e,W,Z,w,z) = ^ (W,Z,w,z) + e 1/2 i/j 1/2 (W, Z, w, z) + e l ^(W, Z, w, z) + (2.9) 

with 

S(e) = So + e 1/2 ^i/ 2 + e 1 ^ + (2.10) 

We substitute these expressions into ()2.7|) and solve the resulting equation order by order in 
powers of e 1 / 2 . 

Note: The description in this section is purely formal. In particular, it does not take into 
account the cutoffs that are necessary for rigorous results. The mathematical details are 
dealt with in the next section. 

Order The order e° terms require 

[h(e, W, Z) - E(e, W, Z)} </>o + E ^ = S ^o- 
We solve this by choosing 

So = Eq, 

and 

MW, Z, w, z) = f (W, Z, w, z) $(W, Z), 

where <&{W, Z, ■) is a normalized ground state eigenvector of h(e, W, Z). Under our as- 
sumptions, we can choose $(W, Z, ■ ) to be real, smooth in (W, Z), and independent of e. 
This choice satisfies 

( $(W, Z,-), Vw,z*(W, Z, ■ ) ) Hh = 0, (2.11) 

where the inner product is in the electronic Hilbert space. We assume that fo(W, Z, w, z) 
is not identically zero. 

Order 1/2 The order e 1//2 terms require 

[h(e, W, Z) - E(e, W, Z)\ tfj 1/2 + E ^ 1/2 = S Vi/a + £1/2 ^o- 

The components of this equation in the <&(W, Z) direction in the electronic Hilbert space 
require 

£1/2 = 0. 
9 



The components of the equation orthogonal to §(W, Z) in the electronic Hilbert space require 

[h(e, W, Z) - E(e, W, Z) ] V1/2 = 0, 



so 

ij 1/2 (W, Z, w, z) = f 1/2 (W, Z, w, z) ®(W, Z). 
Orders 1 and 3/2 By similar calculations, the order e 1 and e 3//2 terms yield 

£\ = £3/2 = 0, 

V>i(W, Z, w, z) = h(W, Z, w, z) $(W, Z), and 
if> 3/2 (W, Z, w, z) = f 3/2 (W, Z, w, z) $(W, Z). 

Order 2 The order e 2 terms that are multiples of <&(W, Z) in the electronic Hilbert space 
require 

" \ Z ' W ' Z) ~\ 1^ {W ' Z > W ' Z) + Enf{W ' Z) f ° {W ' Z ' W ' Z) 

= £ 2 f (W,Z,w,z), (2.12) 

where E^(w, z) is given by (|2.fij) . 

Because of the form of E^(w, z), (|2.12jl does not separate into two ODE's. We do 
not know £ 2 or / exactly, although accurate numerical approximations can be found easily. 
These eigenvalues and eigenf unctions describe the coupled anharmonic vibrational motion 
of all three nuclei in the molecule. As we commented earlier, hypotheses (|2.3jl or (|2.4jl 
guarantee that the eigenvalues £ 2 are discrete and bounded below, with normalized bound 
states fo(W, Z, w, z) in (w, z) for any (W, Z). 

Later in the expansion, we choose the operator T 3 so that f has no (W, Z) dependence. 
With this in mind, equation (|2.12j) determines £ 2 and a normalized function fo(w, z) (up to 
a phase) for any given vibrational level. 

The terms of order 2 that are orthogonal to Q(W, Z) require 

[h(e, W, Z) - E(e, W, Z)\^ 2 = 0, 

Thus, 

^ = f 2 (W, Z, w, z) $(W, Z). 
10 



We split the scalar functions f a (W, Z, w, z) with a > into two contributions 

f a (W, Z, w, z) = fl(W, Z, w, z) + fu(W, Z, w, z) 

where for each fixed W and Z, fa(W, Z, • ) is a multiple of /o(-, • ), and ft^Wi Z,-,-) 
perpendicular to /o(-,-) in L 2 (M. 2 , dwdz). Furthermore, we choose the operators T 3+m / 2 
later in the expansion so that fa(W, Z, ■ ) has no (W, Z) dependence. We will not precisely 
normalize our approximate eigenf unctions, so we henceforth assume /1(W, Z, w, z) = for 
all a > 0. 

Order m/2 with m > 4 We equate the terms of order m/2 and then separately examine 
the projections of the resulting equation into the <3>(W, Z) direction in the electron Hilbert 
space and into the direction perpendicular to &(W, Z). 

From the terms in the $>(W, Z) direction, we obtain the value of £ m /2 and an expression 
for f( m ~4:)/2(W, Z, w, z) = /(m_4)/ 2 (^ w i z )- When m = 6 we choose T 3 so that f can 
be chosen independent of (W, Z). When m > 6, we choose T m / 2 , so that ff m _ 6 \/ 2 can ^ e 
taken to be zero. 

The terms orthogonal to <&(W, Z) in the electronic Hilbert space give rise to an equation 
for [h(e, W, Z) — E(e, W, Z) } ip m /2- This equation has a solution of the form 

4> m/2 (W, Z, w, z) = (fl /2 {W, Z, w, z) + f± /2 (W, Z, w, z) ) $(W, Z) 

+ ^i /2 (W,Z,w,z), 

where V'm/2 ^ s obtained by applying the reduced resolvent operator [ h(e, W, Z) — E(e, W, Z) ]~ 
to the right hand side of the equation. 

In the next section, we prove that this procedure yields a quasimode whose approximate 
eigenvalue and eigenvector each have asymptotic expansions to all orders in e 1 / 2 . 

3 Mathematical Considerations 

In this section we present a mathematically rigorous version of the expansion of Section |21 
This involves inserting cutoffs and proving that many technical conditions are satisfied at 
each order of the expansion. 
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E T (w, z) 



Proposition 3.1 Assume \2. C J\) or \2.J$ . 

Id 2 Id 2 

Then, the spectrum of Hnf = — - tt~^ ~ t: tt^ + Enf(w,z) is purely discrete. 

2 aw 1 2 oz 2 

Proof We use Persson's Theorem (see, e.g., to show that the essential spectrum of 
H NF is empty. This theorem says that if V G L p (IR n ) + L°°(IR n ), with p = 2 if n < 3, p>2 
if n — 4, and p > n/2 if n > 5, then the bottom of the essential spectrum of H = —A + V 
is characterized by the behavior of the operator at infinity. More precisely, 

inf a ess (H) = sup inf ( ^ H J ] : <p G C^(M n \ K) 
k compact 

Since E^p does not satisfy the hypotheses, we replace it with a cut off potential 

E NF (w, z), if E NF (w, z) < T, 
T, otherwise. 

The operator Hp = — A + E F is self-adjoint on the domain of — A. Because is a core for 
both Hnf and Hp, and ( (p, Hnf V 9 ) > ( V 9 ; Ht V? ), for any ip G C£°, the min-max principle 
shows that 

inf a eS s(H NF ) > inf a ess (H T ). (3.1) 

Under hypothesis (|2.3|) or ()2.4|) with a 2 > 0, E^f is arbitrarily large for all large argu- 
ments. Persson's Theorem easily shows that inf <7 e ss(-f^r) = T for all large positive T. 
Inequality ()3.1|) immediately implies the proposition. 

Thus, we need only consider the case of hypothesis ()2.4|) with a 2 = 0. Since Et < T, we 
see that inf cr ess (HT) < T. We shall use Persson's Theorem to prove the reverse inequality. 

Consider a square -ft'(-R) of side 2R > 0, centered at the origin, and let p> G C^(M n \ 
K{R)). We observe that E NF (w, z) = T on the set L R = {w < -R} U {\z\ > R, \w\ < R}, 
for R > 2ai/ci3, provided T < a\R 2 . Therefore, 

/ <p(w,z) (-dl/2-d 2 z /2 + E T {w,z)) p{w,z)dwdz > / T \(f{w, z)\ 2 dw dz. 

JL R J Lb. 

For w > R, we estimate the integral as follows 

/ <p{w,z) {-dl/2 - <9 2 2 /2 + E T {w, z)) (p(w, z) dw dz 

J{w>R, zeR} 



> / / (p(w,z) (—dl/2 + E T (w,z)^ (p(w,z)dzdw 

u>R JzeR 



((f{w, •), hr(w)(f(w, -)} 2 dw. (3.2) 
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For each value of w > R, the operator Iit{w) is a one-dimensional Schrodinger operator with 
potential given by a (cut off) symmetric quartic double well. Hence hxiw) always has a 
ground state below T, for any large T. We shall show that the ground state converges to T 
as w — > oo. 

To do this, we show that hj<(w) — > — <9 2 /2 + T in the norm resolvent sense. By the 
resolvent identity, this follows if we show that \\(Et{w, z) — T)(— <9 2 /2 + tends to zero 

as w — > oo. To show this we use the following claim: 
There exists c(T), such that, 

||V(_ 5 2 + T )-i|| L2(RHL2(R) < C (T) \\V\\ 2 , 

for all V G L 2 (1R). 

By Theorem IX. 28 in [H], given any a > 0, there exists 6 > 0, such that 
IMloo < Q ||( — d^)if\\2 + & 1Mb) for any ip in the domain of —d 2 z . This implies that 

II V {-dz + \\l 2 (R)^L' 2 (R) < || U ||l°°(rhl 2 (r) II (-$2 + |U 2 (R)^L°°( J K) 

< c(T) ||V|| 2 , 

for some finite c(T). This proves the claim. 
In our case, this yields 

\\(E T (w,-)-T)(-d 2 z /2 + Tr 1 \\ < c(T) || (E T (w, ■ ) - T) || 2 , 

where c(T) is independent of w. 

Under the assumption that R > yjT/a\, we have 



(E T (w,-)-T)\\l < [ T 2 dz 

J\w-z 2 a 3 /(2a 1 )\< y /T/a 1 



^2 



2T " (V W + V^M - yfw-y/T/~a 1 



m2 /2T 

V «3 



c(T) — =, as w — > oo. 



Thus by the resolvent formula and Theorem VIII. 19 of JU], hr{w) converges in norm resol- 
vent sense to (-<9 2 /2 + T). 
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We let Pa{H) be the spectral projector on the interval A G M for a self-adjoint operator 
H . Then Theorem VIII. 23 of jTU] implies that for any positive a < b < T, 

lim \\P (aib) (h T (w)) - P (a , b) (-d 2 z /2 + T)\\ = lim \\P {a , b) (h T (w))\\ = 0, 

since a(—d 2 /2 + T) = [T, oo). Therefore, the bottom of the spectrum of ft/r(u>) satisfies 

T > bT{w) = inf a(hr(w)) — > T as 10 — > oo. 

Using this in (|3.2)l . we obtain 

/ <p(w,z) (~dl/2 - d 2 /2 + E T (w, z)) <p(w, z) dw dz 

J{w>R, z£i?} 



> inf bx{v) / |(/3(w,z)| 2 ciwtiz. 

V ^ R J\w>R,zGR\ 



'{w>R,z£R} 

Combining all the estimates, we see that for any ip G C^(IR n \ K(R)) 



(ip, H T ip) > inf bxiy) || ip 

v>R 



|2 



provided ai R 2 > T. Thus, 



inf aess(Hr) > sup inf 6r(t>) = T. 



Since T is arbitrarily large, this implies the proposition. I 

In the usual Born-Oppenheimer approximation, the semiclassical expansion for the nuclei 
is based on Harmonic oscillator eigenfunctions. They have many well-known properties. 
Our expansion relies on the analogous properties for eigenfunctions of H^p. The following 
proposition establishes some of the properties we need in an even more general setting. 

Proposition 3.2 Let V be a non-negative polynomial, such that H = —A + V has purely 
discrete spectrum. Let <p(x) be an eigenvector of H , i.e., an L 2 (R") solution of Hip = Sip, 
where £ > 0. Then, ip G C°°(R n ) and Vip G L 2 {lR n ). Moreover, for any a > 0, 

peD{e a{x) ), Vip G D{e a{x) ), and Aip G D{e a{x) ), 

where (x) 



1 + x 2 , and D(e a ^) denotes the domain of multiplication by e a ^ . 
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Proof Since V G C°°, elliptic regularity arguments (see e.g., [7], Thm 7.4.1) show that all 
eigenf unctions are C°°. 

We first show that the W<p is L 2 . Since V > 0, the quadratic form defined by 

on = Q( — A) fl Q(V), is closed and positive. Here Q(A) means the quadratic form 

domain of the operator A. Since D(H) C Q(h), any eigenvector of H belongs to 

Q(-A) = { G L 2 (lR n ) : ||V^|| < oo }. 

Thus, Vy? G L 2 . 

Next, we prove y? G D(e a( - X ^), for any a > by a Combes-Thomas argument, as presented 
in Theorem XII. 39 of ^2]- We describe the details for completeness. Let a G M, and let v 
denote Xj for any j G {1, • • • , n}. We consider the unitary group W(a) = e iav for a <E M, 
and compute 

H(a) = W(a) (-A + V) W(a)- 1 = H + iad v + a 2 . 

The operator id v is i?-bounded, with arbitrary small relative bound, since V > 0. Thus 
{if (a)} extends a self-adjoint, entire analytic family of type A, defined on D(H). We note 
that since H(0) = H has purely discrete spectrum, its resolvent, -Ro(^) is compact, for any 
A G p{H) = (D \ <j(H). Hence, R a {^) = {H(a) — A) -1 is compact for any a G M, and 
hence, for all a G (D, if A G p(H(a)). It is jointly analytic in a and A. The eigenvalues 
of H(a) are thus analytic in a, except at crossing points, where they may have algebraic 
singularities. Since for a real, W(a) is unitary, the eigenvalues are actually independent of 
a, and a(H(a)) = &(H), for any a. 

Let P be the finite rank spectral projector corresponding to an eigenvalue S of Hnf • 
Then, for a G M, P(a) = W / (a)PW / (a)" 1 is the spectral projector corresponding to the 
eigenvalue S of H(a). By Riesz's formula and the properties of the resolvent, P(a) extends 
to an entire analytic function that satisfies 

W(a )P(a)W(a )- 1 = P(a + a). 

for any cto G M. 

By O'Connor's Lemma (Sect. XIII. 11 of ^2])> this yields information about the eigenvec- 
tors. If if = Pip, the vector (p a = W(a)(p, defined for a G M has an analytic extension to the 
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whole complex plane, and is an analytic vector for the operator v. Therefore, cp G D(e a ^), 
for any a > 0. By taking all possible x,'s for v, and noting that D(e a ^) = D(e a ^i \ x ^"j } we 
see that (p G D(e a{x) ). 

From this, it follows that Ay? G D(e a ^) for any a > as well, since for any 5 > 0, 



/ 




|Av,(x)| 2 






s 


e 2a(x> 


1 Wx) - 


-S)(p(x) | 2 




II (v 




e" 50 II 

° || oo 


|| e («+5/2){-> 


*>(•) ir 


oo. 











Finally, Lemma 13.31 below shows that Vy? G D(e a ^). To apply this Lemma in our 
situation, we let p(x) = e a ^ and note that for any a > 0, 

(Ve a(x,) )/e a(x) = aV(x) = ax/(x) 

is uniformly bounded. I 

Lemma f3 .31 requires some notation. Letting p(x) be a positive weight function, we intro- 
duce the space 

F 2 W = if : \\f\\ 3 n = / (|/(x)| 2 + |A/(x)| 2 ) < oo). 

We write ||/|| 2 = J Rn \ f(x)\ 2 p(x) dx, for any / G L 2 (W l ,p(x)dx), and ||/|| 2 = J Rn |/(x)| 2 rfx 
when the weight is one. 

Lemma 3.3 Let p G C 1 be positive, and assume that there exists a constant C < oo, such 
that \(Vp(x))/p(x)\ < 2C for all x G R n . Then, for any f G F 2 

||V/|U < C\\f\\ w + VII/IU P/C + (3.3) 



We present the proof of this technical lemma in Section |3J 

We now state and prove the following Corollary to Proposition 13.21 
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Corollary 3.4 Assume the hypotheses of Proposition \S.2l Let R(X) be the resolvent of 
H = —A + V for A ^ o~(H), and let Pg be the finite dimensional spectral projector of H on 
£. Let r(S) = ( {H — £)\(i-p s )i? ) 1 be the reduced resolvent at £. Then, e a ^ R(X) e~ a ^ 
and e a ^ r(£) e~ a ^ are bounded on L 2 (lR n ). 

Proof We use the notation of the proof of Proposition ^. 21 We know that R a {\) is compact 
and analytic in a G C, if A ^ cr(H). Hence, for any tpi, ip 2 G Co°, the map from M x p(H) 
to (C given by 



is uniformly bounded by C [| Vx|| II^H on any given compact set of M x p(H) for some C. 
From this we infer that for any a > 0, e a ^ R{\) e~ a ^ is bounded in L 2 (]R n ), uniformly for 
A in compact sets of p(H). Since the reduced resolvent r{£) can be represented as 



where Cg is a loop in the resolvent set encircling only £, the boundedness of e a ^ r(£ ) e a ^ 



To show that the terms of our formal expansion all belong to L 2 , we use the following 
generalization of Proposition 13.21 We present its proof in Section 

Proposition 3.5 Assume the hypotheses of Proposition and let if be an L 2 solution of 
(—A + V — £) cp = 0. Then, for any a > 0, and any multi-index a e N n , D a ip 6 D(e a( - X ' ) ), 
where D a = d" 1 d? 2 ■■■ d? n . 

We now prove that our formal expansion leads to rigorous quasimodes for the Hamilto- 
nian H\{e) given by (|2.1|) . Theorem 13.61 summarizes this result for the leading order, while 
Theorem 13. 71 handles the arbitrary order results. 

Theorem 3.6 Let h(e, W, Z) be defined as in Section\^with W shifted so that Wo = 0. We 
assume h(e, W, Z) on 7i e i is C 2 in the strong resolvent sense for (W, Z) near the origin. 
We assume its non- degenerate ground state is given by 



(a, A) 



(^,e av R (\)e- av ij 2 ) 




(3.4) 



follows. 




e, W, Z) + S{e, W, Z), 



(3.5) 
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under hypothesis \2.'J\) or \2.J$ , and we denote the corresponding normalized eigenstate by 
§(W, Z). Suppose the remainder term S is uniformly bounded below by some r > — oo and 
that \S\ satisfies a bound of the form 

| S(e, W, Z)\ < C \W a Z 2f3 \ (3.6) 

a+/3>3 

for (W, Z) in a neighborhood of the origin. Here C is independent of e, the sum is finite, and 
a and (3 are non-negative integers. Let fo{w,z) be a normalized non-degenerate eigenvector 
of H NF , i.e., 

(-d 2 w /2 - d 2 z /2 + E NF (w, z)) f = S 2 f , 

with 

E NF (w,z) = a\ w 2 + (^a 2 — a 3 w^j z 2 + a 4 z 4 . 
Then, for small enough e, there exists an eigenvalue £{e) of H^e) which satisfies 

S{e) = E Q + e 2 S 2 + 0(e«), 

for some £ > 2 as e — > 0. 



Remarks 1. At this level of approximation, it is not necessary to require the eigenvector 
$ to satisfy condition (j2.11j) or to require h(e, W, Z) be real symmetric. 
2. We have stated our results for the electronic ground state, but the analogous results 
would be true for any non-degenerate state that had the same type of dependence on e. 

Proof: In the course of the proof, we denote all generic non-negative constants by the same 
symbol c. 

Our candidate for the construction of a quasimode is 

* (e, W, Z) = FiW/e^) F{Z/e 5 *) f (W/e, Z/y/e) $(W, Z), (3.7) 

where F : M — > [0, 1] is a smooth, even cutoff function supported on [—2, 2] which is equal to 
1 on [—1, 1]. One should expect the introduction of these cutoffs not to affect the expansion 
at any finite order because the eigenvectors of H\(e) are localized near the minimum of 
E\(e, W, Z). Thus, the properties of the electronic Hamiltonian for large values (W, Z) 
should not matter. The choice of a different cutoff for each variable is required because 
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these variables have different scalings in e. We determine the precise values of the positive 
exponents Si and 5 2 in the course of the proof. We also use the notation 

F(e,W,Z) = F(W/e Sl ) F(Z/e 52 ). (3.8) 

We first estimate the norm of ^/q. 



\*q\\ 2 = [ \F(e,W,Z) f (W/e, Z / V~e)\ 2 \MW, Z)\\ 
Jr 2 

\fo(W/e, Z/^e)\ 2 dWdZ 



[ (l-F 2 (e,W,Z)) \f (W/e, Z/^)\ 2 dWdZ. 
Jr 2 



'R? 

The first term of the last expression equals e 3 / 2 , by scaling, since / is normalized. If Si < 1 
and S 2 < 1/2, the negative of the second term is bounded above by 



L>A \fo(W/e,Z/V~e)\ 2 dWdZ 

J R~ 



\Z\>e S 2 

= e 3 / 2 f w>el _ Sl e- 2a ^+N) e MW+\*\) | /o ( W) z )f dwdz 

\z\>S/ 2 ~ S 2 

< e 3/2 e -2o(lA( 1 -*i)+l/ £ ( 1 -*2)) || e a(H+H)y o ||2 

= 0(0. 

since / G D(e a <W z )>). Hence, 

||^ Q || = e 3/4 (1 + 0(e°°)), where the 0(e°°) correction is non-positive. (3.9) 
Next we compute 

(Hi(e) - (E + e 2 S 2 )) * g (e, W, Z) (3.10) 
= S(e, Z) /o^zJI^^f.Z) $(W,Z) 

Y^ + ^l) .F(e,W,Z)$(W,Z)) / (w,*)kz, 
-e 3 (0„,/ o (u;,*)kz) d w (F(e, W, Z)$(W, Z)) 
-e>' 2 (dj (w,z)\ w , z ) d z (F(e,W,Z)<S>(W,Z)), 

19 



where we have introduced the shorthand fo(w, z)\w,z = fo{W/e, Z / yfe) and used the iden- 
tity 

E{e, W, Z) - e 2 E NF {W/e, Z/y/e) = 0. 



Also 



1 



d w f(e } W } Z) = ^- F'(W/e Sl ) F(Z/e 52 ' 



€ 



Si 



d z F(e,W,Z) = F{W/e 5 ') F'(Z/e 6 *) 

and, by assumption, ||9^9^$(W, ^)||w el is continuous and of order e° in a neigborhood of 
the origin, for jj, + v < 2. Therefore, 



sup 

R 2 


d w (F(e,W,Z)$(W,Z)) 


llw e i 


< 


c 


sup 

R 2 


d z (F(e,W,Z)$(W,Z))\ 




< 


c 


sup 

R 2 


d? w {F(e,W,Z)$(W,Z)) 


llw e l 


< 


c 


sup 

R 2 


d?(F{e,W,Z)$(W,Z))\ 


\n e i 


< 


c 


where all vectors 


are supported in { (W, Z) 


■ \w\ 


< 





(3.11) 



appears in (|3.10j) . multiplied by one of the scalar functions fo(w, z)\w,z, (d w f (w,z)) \w,z, 
or (d z f (w, z)) \ w,z- In turn, each of these functions belongs to L 2 (IR 2 ) by Proposition 13.21 
and each one has norm of order e 3//4 because of scaling, e.g., 

( I (d w f )(W/e, Z/V~e) ? dW dZ^j = e 3 / 4 || d w f \\ L , {R2) . 

Therefore, the norms of the last three vectors in ()3.1()|) are of order e 3 / 4 times the corre- 
sponding power of e stemming from (|3.11|) . 

We now estimate the norm of the term that arises from the error term S. From our 
hypothesis on the behavior of S, we have 

II S T / $|| 2 = J ]w ^ /eh \f (W/e,Z/V~e)S(W,Z)\ 2 dWdZ 

Z|<2/e <5 2 

E L< 2/ ^ \fo(W/e,Z/Vi)\* \ W°ZV\ 2 dWdZ 



< C 

* ' l\W\<2/ 

a+/3>3 J \Z\<2/e S 2 
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< c J2 e 2 ^ +2 ^ [ \f (W/e, Z/^e)\ 2 dWdZ 



«+/3>3 



f . . e 2(a5 1 +2/35 2 ) £ 3/2 



£ 

a+/3>3 

where the sums are finite. 

Collecting these estimates and inserting the allowed values of a and (3, we obtain 

|| (H 1 {e)-(E + e 2 £ 2 )) * Q || 

< c e 3/4 ( e 351 + e 2(<5l+<52) + e Sl+A52 + e &&2 + e 4 ~ 2Sl + e 3 ~ 2S2 + e 3 ~ 51 + e 5/2 ~ 5 ' 2 ) 

We further note that 8 X < 1 and 5 2 < 1/2 imply e 4_2<51 < g 3- * and e 3 ~ 252 < e 5 / 2 " 52 . This, 
together with (|3.9|) . shows that for small enough e, 

|| ~ {Eg + e 2 ^)) || < c / e 35! , e 2(5!+5 2 ) , e 5i+45 2 , e 6«5 2 , £ 3-5i , e 5/2-5 2 \ 

II *Q I " ^ ' 

We still must show that all terms in the parenthesis above can be made asymptotically 
smaller than e 2 . This can be done if there exist choices of 5\ and 5 2 such that all exponents 
in the parenthesis above are strictly larger than 2. The inequalities to be satisfied are 

< Sx < 1, Sx> 2/3, 5t + 5 2 > 1, 0<5 2 < 1/2, 5 2 > 1/3, <5i + A5 2 > 2. 

Satisfying these is equivalent to satisfying 

2/3 < 5t < 1 
1/3 <5 2 < 1/2 

which defines the set of allowed values. The best value, 

f = max min {35 u 2{8 1 + 8 2 ), 5i + 45 2 , 65 2 , 5/2 - 5 2 } > 0, 

0<« 1 <l 
0<« 2 <l/2 

is obtained by straighforward optimization and is given by £ = 15/7, obtained for 5/7 < 
<5i < 6/7 and <5 2 = 5/14. With such a choice, there exists an eigenvalue £(e) of H\{e) that 
satisfies 

£(e) = E + e 2 £ 2 (e) + 0(e*), 
with £ = 2 + 1/7. ■ 

We now turn to the construction of a complete asymptotic expansion for the energy level 
£(e) of #i(e), as e -> 0. 
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Theorem 3.7 Assume the hypotheses of Theorem Iff. 61 with the additional condition that 
h(e, W, Z) on 7i e i is C°° in the strong resolvent sense in the variables (e, W, Z). Then the 
energy level £(e) of -ffi(e) admits a complete asymptotic expansion in powers of e l l 2 . The 
same conclusion is true for the corresponding quasimode eigenvector. 

Proof Our candidate for the quasimode is again the formal expansion ([2.9)1 truncated at 
order e N l 2 and multiplied by the cutoff function ([3.8)1 . i.e., 

N 

* Q (e, W, Z) = Tie, W, Z) ^ e^ 2 ^ /2 (W, Z, W/e, Z/V^)- 

i=o 

We shall determine ?/>j/2 and Tj/2 in (|2.8jl explicitly, but first we introduce some notation 
for certain Taylor series. Expanding in powers of e 1 / 2 , we write 

T j/2 (W, Z) - T j/2 (ew, t^z) 
= T j/2 (W, Z) - T j/2 (0, 0) - e 1 / 2 9 z T i/2 (0, 0) + e (d w T j/2 (0, 0)w + d z ,T j/2 (0, 0)w 2 /2) + ■■■ 

oo 

= T j/2 (W, Z) - T j/2 (0, 0) + J2 *) ^- 

k=l 

Next, our hypotheses imply that the function S(e, W, Z) in ([3.5)1 

Ex(e,W,Z) = E + E(e, W,,Z) + S(e, W, Z) 
is C°° in (e, W, Z). Using flUD), we write 

oo 

Ex(e, ew, e 1/2 z) = E + e 2 E NF (w, z) + ^ e m ' 2 S m/2 {w, z). 

m>6 

Note Because we have assumed Ei(e, W, Z) is even in Z, S m / 2 (w, z) = when m is odd, 
but the notation is somewhat simpler if we include these terms. 

We use this notation and substitute the formal series (|2.9|) and (j2.1(J|) into the eigenvalue 
equation ([2.7)1 . with H 2 given by ([2.8)1 . For orders n/2 with n < 4, we find exactly what we 
obtained in Section |21 When n > 5, we have to solve 

[h(e,W,Z)-E 1 (e,W,Z)]^ n/2 (3.12) 

+ E NF (W, z)^(n-4)/2 + S 6/2 (W, 2;)-0(„- 6 )/2 + SV/ 2 (w, ^)^(n-7)/2 H h S n / 2 (w, z)^o 
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+ (Tn (W, Z) - Tn (0, 0))^0 - rffi (W, Z)^ - (W, z)^ T„ 



( "~ 6 ) 

2 (w, z)ip< 







+ (T^(W, Z) - IWO, 0))^i/2 - r^_% /2 (w, z)^j 



T, 



w,z)if)y 2 



+ (TV (W, Z) - TV (0, 0))fe - («,, z)^ {n _ 7)/2 

2 2 2 2 

+ (T f (^,Z)-T|(0,0))Vv-6)/2 

= £2 1p(n-4)/2 + £5/2 '0(n-5)/2 + " " " + ^n/2^0, 

with the understanding that the quantities 5*, T and r that appear with indices lower than 
those allowed in their definitions are equal to zero. 
We solve (j3.12l bv induction on n. We assume that 

%2, ipj- /2 {W,Z,w,z), T j/2 {W,Z) for j<n-l, and 
fj/ 2 {W, Z, w, z) for j < n - 5 

have already been determined, with fh 2 (W, Z, w, z) = 0, for j > 1. 

We project (j3.12|) into the <&(W, Z) direction and the orthogonal direction in the elec- 
tronic Hilbert space to obtain two equations that must each be solved. 

First, we take the scalar product of (j3.12|) with $>(W, Z) in the electronic Hilbert space 
to obtain 

E NF (w,z)f( n - 4 )/2 + S 6 / 2 {w,z)f( n -6)/2 + S 7 / 2 (w, z)/( n -7)/2 + + S n / 2 (w,z)f 

+ (Tn(W, Z) — T"(0, 0))/ - S\{w,z)h - S 2 (w,z)f - ••■ - r^\w,z)f 

2 * 2 2 2 

+ (T^(W, Z)-T v (0,0))/ 1/2 ^ T^l )/2 {w,z)f l/2 - - rf^t^/i/a 



+ (T ; (WiZ)-r ; (0,0))/»-r - ^Kz)/ M/J 



+ {Te(W, Z) — Te(0, 0)) /( n _ 



6)/2 
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77 ^w,zf(n-i)/2 — 



- \ ($(W, Z), d\ z ^ m ) Uel - l - (<S>(W, Z\ d 2 mw ^ 8)/2 ) Hel 

— £2 /(n-4)/2 + ^5/2 /(n-5)/2 + ' ' " + ^n/2 fo- (3.13) 

We further project ()3.13j) into the fo direction and the orthogonal direction in L 2 (1R 2 ) to 
obtain two equations that must each be solved. 

We take the scalar product of (|3.13|) with f in L 2 (R 2 ). Using /jL = for j > 1 and 
(-§A^ 2 + E NF (w, z)) f = £ 2 f , we obtain 

£ n /2 = T f (W,Z)-T f (0,0) (3.14) 

ra n— 7 n— 6— fe 

+ (Jo, 5^72 /(j-6)/2)i2( R 2) - S (/o, T j y 2 " 0+fc))/2) /fc/2)L2(R2) 

j=6 fc=0 j=l 



fo, [{*{W, Z), d% z ^ 5)/2 ) Hel + (*{W,Z), d 2 w ^ in . 6)/2 ) Hel 

+ 1 (*(W,Z), 9l^ (n _ 6)/2 ) Wei + i (SflV.Z), ^,^(n-8)/2> WeIJ , 

We can solve this equation for S n / 2 if the right hand side is independent of (W, Z). This will 
be true, if we choose 

n Ti—7 n— 6— k 

Tn(W, Z) = / {j - 6 )/2>l W + ^ £ </o, r ( ;- (j+fe))/2) / fc/2 ) L2(M2) 



+ (/o, ($(w, z), dl^ {n . 5)/2 ) Hel + mw, z), d 2 w ^ {n ^ )/2 ) 



We then are forced to take 

£ n/2 = -T™(0, 0). 

The first non-zero Tj/ 2 (W, Z) is 

T 6/2 (W, Z) = l - ($(W, Z), d 2 z ^W, Z)) Hel - (fo, S 3 fo)m^y (3.15) 

So, 

£3 = (fo,S 3 f ) L 2 im - \ ($(0, 0), (<9f$)(0, 0)) WeJ . 
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We next equate the components on the two sides of (|3.13jl that are orthogonal to /o in 
L 2 (IR 2 ). The resulting equation can be solved by applying the reduced resolvent r NF (£ 2 ), 
which is the inverse of the restriction of (— |A Wj2 + E^p — £2) to the subspace orthogonal to 
fo. We thus obtain 



/(n-4)/2 — r NF (£ 2 



E £ (ri-j)/2 fj/2 - E ( S j/z( W ' Z ) f(i-e)/2) ± 

- i=i i=6 

n— 7 n— 6— k n— 6 

+ E E (^^(^^/va) 1 + E^-i)/^ ' )-^)/^^^))/^ 

fc=o i=i i=i 
+ ( Z), d\ z * (n _ 5)/2 )^ ei + ( Z), fl^ * (n _ 6)/2 >^ 



+ ( $(W, Z), * (n _ 6)/2 )^ + ( Z), ^ * (n _ 8)/2 )^ ei 



This solution has /(„_4)/ 2 = fL_±\/ 2 orthogonal to / , as claimed in Section [21 The first 
non-trivial fj/2, for j > 1 is 



ft{W, Z, z) 



r N F{£2) (S 3 (w,z) f {w,z)) L . 



(3.16) 



Next, we equate the components of (|3.12|) that are orthogonal to <&(W, Z) in TC e [. We 
solve the resulting equation for T/>„/ 2 by applying the reduced resolvent r(W, Z) of /i(e, W, Z) 
at £ x (e, W,Z). This yields 

n— 5 n 

r(W,Z) ^ S {n -j )/2 ipj- /2 - S j/^w,z) Vy_ 6 )/2 ( 3 - 17 ) 



/2 



i=i 

n— 7 n— 6— fc 

E E 

fe=o i=i 



i=6 



n— 6 



(n-(j+k))/2)/ 
T i/2 ^ 



, *)^/ 2 + £(T (n _ j)/2 (0, 0) - T (n _ i)/2 (W, Z))^ /2 
j'=i 



+ i d W,W AnS)^) 1 ~ {^ A w,z + E NF (W,Z) - £ 2 ) ^_4)/2 

The first non-zero component ip^ 2 with j > 0, is 

^ /2 (W, Z, w, z) = (dj )(w,z) r(W, Z) (d z <f>)(W, Z). 



(3.18) 



Finally, Proposition l3.9l below shows that each ipj/ 2 in this expansion belongs to D(e a ^ w " e+ \ z \' v ^- ) ). 
As a result, whenever a derivative acts on the cutoff, it yields a contribution whose L 2 norm 
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is exponentially small. This way, we can neglect such terms. For example 

(dwFie, W, Z))ij j/2 (W, Z, W/e, Z/y/e) = e^ 1 F\W/^) F(Z/&) ^ /2 {W, Z, W/e, Z/^e). 

The square of the L 2 norm of this term is bounded by a constant times 

^^(^^W/e^/V^I^^^I^+l^/^e-^ 1 -' 1 )^^ = 0(e°°). 



'\W\/e d i>l,\Z\/e"2<2 

Thus, we have constructed the non-zero quasimode ()3.7|) that satisfies the eigenvalue 
equation up to an arbitrary high power of e 1 / 2 . I 

The proof of Proposition 13.91 relies on the following lemma. 

Lemma 3.8 Let V be a polynomial that is bounded below, such that the spectrum of H = 
-|A + V purely discrete. Let tp G C°°(IR n ) satisfy D a tp G D{e a ^), for all a G N n and any 
a > 0. If R(X) denotes the resolvent of H , then D a R(X)(p e D(e a ^) for all a G N n and all 
A in p(H). The same is true for D a r(£)ip, where r(S) is the reduced resolvent at S. 

Proof We first note that elliptic regularity implies that the resolvent -R(A) maps C°° func- 
tions to C°° functions. Next, applied to smooth functions in L 2 , we have the identity 

[d Xj ,R(X)} = R(X) (d Xj V)R(X). 

We claim that the operators on the two sides of this equation have bounded extensions to all 
of L 2 . To see this, note that D@V is relatively bounded with respect to V for any (3 G N n , 
because V is a polynomial. Furthermore, since H > V, we see that D@V is relatively 
bounded with respect to H, which implies the claim. Hence, for tp as in the lemma, we have 

d x . R(X) il> = R(X) d x . <p + R(X) (d x .V) R(X) ip. (3.19) 

The first term on the right hand side of this equation belongs to D(e a( - X ^) since R(X) maps 
exponentially decaying functions to exponentially decaying functions (see Corollarv l3.4|) . The 
same is true for the second term, with a possible arbitrarily small loss on the exponential 
decay rate, due to the polynomial growth of d x .V. This provides the starting point for an 
induction on the order of the derivative that appears in the conclusion of the lemma. 

We now assume that for some a G N n , D a R(z)(p is a linear combination of smooth 
functions of the form R(z)(D^V)R(z) ■ ■ ■ R{z)(D^- 1 V)R(z)D^^ all of which belong to 
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D(e a ^), for any a > 0. We assume every 7 that occurs here has I7I < |a|. Then d Xj D a R(z)(p 
is a linear combination of elements of the form 

d x .(R(z)(D^V)R(z) ■ ■ ■ R(z)(D^V)R(z)D^ ( p) 

Applying ()3.19j) successively, we see that the structure is preserved. Since all D lk V are 
polynomial, Corollary 13 . 41 implies the result. 

The statement for the reduced resolvent follows from the representation ()3.4j) . I 



Proposition 3.9 Assume the hyptheses of TheoremUT} Leti/jj/2(W, Z, w, z) be determined 



by the construction above, where (W, Z) belongs to a closed neighborhood Q of the origin and 
(w, z) G 1R 2 . Then ipj/ 2 is C°° , and the function G(w, z) = sup^^g^ \ipj/ 2 {W, Z,w, z)\ 
belongs to D(e a{{w > z)) ). 

Proof The hypothesis on the Hamiltonian and the properties of the normal form Hnf 
proven above imply that Z) and r(W, Z) are smooth, and that r^i?^) maps smooth 

functions to smooth functions. We also know that the non-degenerate eigenstate fo is smooth 
and belongs to D(e a ^ w,z ^). The smoothness of ^^(W, Z, w, z) follows trivially in VL x IR 2 . 

Concerning the exponential decay, we observe that the (w, z) dependence of ipj/2 stems 
from the successive actions of derivatives, reduced resolvents, and multiplications by polyno- 
mials in (w, z), acting on the eigenstate fo. Lemma f3.8l applied in conjunction with Propo- 
sition 13.51 shows that the exponential decay properties are preserved under such operations. 



4 Technicalities 

In this section, we present the proofs of Lemma 13.31 and Proposition 13.51 

Proof of Lemma 13.31 We first note that the hypothesis on p implies p(x) > for any 
16K", and that 

Let B R GR n bea ball of radius R > 0. We first show that / G L 2 (B R+1 ) and A/ G L 2 (B R+1 ) 
imply / G H 2 (B R ), where 

H 2 {B R ) = { f G L 2 (B R ), Vf G L 2 (B R ), and Af G L 2 (B R ) }. 
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We denote the usual H 2 (Br) norm by || • \\h 2 (b r )- We now show the existence of a constant 
K(R) > 0, which depends only on R, such that 

/ |V/| 2 < K(R) I (|A/| 2 + |/| 2 ). (4.2) 

JB R JB R+1 

Note This estimate does not hold in general if the balls over which one integrates have the 
same radius. 

We set g = Af on B R+ i and g(x) — if > R+ 1. We can then decompose / = f\ + f 2 
with fi and f 2 solutions to 

Afi = g, fi\dB R+3 = 
Af 2 = \x\<R + l. 

Thus, /i G H 2 (B R+3 ), and there exists a constant Ci(R), which depends only on R, such 
that 

\\fi\\ H *(B R+3 ) < ci(i2) \\AM\ L 2 (Br+3) = Cl (R) \\Af\\ LHBR+l) , (4.3) 
so that 

l|V/l|| L 2 (Bil+l) < \\fl\\ H 2(B R+3 ) < Cl(i?) \\Af\\ L 2 {BR+l) . 

By the mean value property for harmonic functions, f 2 also satisfies estimate (J4.2)) . for some 
constant K 2 (R) with A/ 2 = (see e.g., Chapter 8 of |Tj). Combining these arguments, we 
see that for c 2 (R) = d(R) + K 2 (R), 

[ |V(/ 1 + / 2 )| 2 < c 2 (R) [ (|AM 2 + |/ 2 | 2 ) < c 2 (R) [ (|A/| 2 + 2(|/| 2 + |A| 2 )). 
Jb r Jb r+1 Jb r+1 

But Ib r+1 l/il 2 - ll/i II H2 (Bfl+3) , so (jOl) implies that (jOl) holds for some constant fC(i?). 

Because of (J4.1j) . we can insert the weight p into this estimate to establish the existence 
of another constant K(R), which depends only on R, such that 

/ p |V/| 2 < K(R) [ p(|A/| 2 + |/| 2 ). (4.4) 
JB R Jb r+1 

In other words, p 1 / 2 V/ £ Lf oc ii f £ F 2 . 

A second step consists in showing that p 1 ^ 2 V/ is in L 2 (IR n ) and satisfies (|3.3|) . Let 
X_r G C°°(]R n ) be a truncation function such that < < 1, with Xr( x ) = 1 if \x\ < R, 
and Xr( x ) =0 if \x\ > R + 1. We can take xr so that || Vxr||oo is independent of R. Let 
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/ G F*, and set f R = xrJ ■ Since V f R = XrVI + f^XR, we see that \\p 1/2 Vf R \\ L 2 {BR) = 

||p 1/2 V/|| L2(Sii) , and 

lim || p V2 (/i? _ / )|| 

H—>oo 

by Lebesgue dominated convergence. By the same argument with A/# = XrA/ + J'Axr + 
2V XR Vf, 

lim |b 1/2 (A/ - (xflAZ + ZAxfl))!!^^) = 0. 

it—* oo 

We have the estimate Ib 1 ^ 2 Vx_r V/||^ 2 ( Bfl+1 ) < c 2 Jb r+1 \b r P 1 ^ 2 I V/| 2 , for some constant 
C2, independent of R. We can cover the set B R+ i\B R by a finite set of balls {i^i^')}^!,...^^), 
of radius 1, centered at points Xj such that \xj\ — R + 1/2. In each of these balls -Bi(j), we 
can apply (|4.4jl (with a constant K\, independent of R), to see that 

/ P |V/| 2 < c* £ / P(|A/| 2 + |/| 2 ), 

JB R+1 \B R j=1 JB2U) 

where B 2 (j) has radius 2 instead of 1. Using LiJ^ B 2 (j) C -Br+3 \ -Br-3, and taking into 
account that certain points are counted (uniformly) finitely many times in the integral, we 
eventually obtain 



Ip^VxrV/HI^j < c 3 f _ p(|A/| 2 + 



where C3 is uniform in R. By the dominated convergence theorem again, this integral goes 
to zero as R goes to infinity. So, we finally obtain 



lim lb 1 / 2 (Af - Af R ))\\ L2{Rn) = 0. 

it— ► OO 

Since f R belongs to H 2 {R n ), the set of compactly supported functions in H 2 , we compute 
V • (v h V/r) = P f R Af R + \Vf R \ 2 P + f R V P V/r. 
Since has compact support, Stokes Theorem and our hypotheses on Vp show that 
P\fn\ 2 = 

< (/ \Af R \* p y (/ |/i ? | 2 P ) 1/ ' + 2C (/ I/rI 2 ^^ (I |V/h|V 



/r Af RP + J f R Vp Vf R 

1/2 / » \ 1/2 / „ \l/2/ r \ 1/2 



or, in other words, 

liv/^ii 2 < ||a/r|u h/bIU + 2c ||v/r|u ll/jilU- 
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This estimate implies f!3.3p for The right hand side of that estimate has a finite limit as 
R — > oo with / in place of on the right hand side. Since 



V |V/| 2 < / p \Vf R \\ 

Br J 

we deduce that p |V/| 2 G L 1 (M n ) and satisfies (Q. ■ 

Proof of Proposition IH751 We use the following Paley-Wiener theorem, Theorem IX. 13 
of HI]: 

Let f G L 2 (M n ). Then e a ^f G L 2 (M n ) for all a < a' if and only if f has an analytic 
continuation to the set {p : \Im p\ < a'} with the property that for each t G M n with \t\ < a' , 
f '(• + it) G L 2 (lR n ) , and for any a < a 1 , sup| 4 | <a / ||/(- + it) || 2 < oo. 

We refer to the conditions on / in this theorem as "the Paley-Wiener conditions." 

Since e a ' x ' (p G L 2 (IR n ) is equivalent to <p G D^e^), Proposition 13.21 shows that <p is 
analytic everywhere and satisfies the Paley-Wiener conditions. The functions p t— > pj (pip) 
and p i— > Pj <pip) also satisfy these conditions. 

As a preliminary remark, we note that for any fixed t G M n , there exist K{t) > Kit) > 
and R(t) > 0, such that if p G M n satisfies ^ ■ p 2 > -R(t), then 



K(t) ^ v: < 



< a-w E p?- ( 4 - 5 ) 

i=i 



i=i 

So, if B R is a ball of radius R with center at the origin, —A<p satisfies 

/ (X>?) \<fip + tt)\ 2 dp <oo, (4.6) 

Jm\B R(t) \ j=1 J 

uniformly for t in compact sets of M n . 

We now start an induction on the length \a\ of the multi- index a in D a ip. We first 
show that p t— > PjPk<fip) satisfies the Paley-Wiener conditions for any j, k G {1, ••-,n}. 
Note that we only need to prove estimates for large values of the |pj|'s. Also, note that if 
Y^j=\Pj ^ R{t) > 1) there exists a constant C(t) > such that 

n 

I (pj + itj) ip k + it k ) | < C(t) Pi ( 4 - 7 ) 

j=0 
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Therefore, ()4.6|) implies that 

/ I (Pj + ttj) {Pk + itk) I 2 \<p(p + it)\ 2 dp 

Jn"\B m 

< c 2 (t) [ dtp*) mp+tt)\ 2 d P 



(t) \ j=l 



< OO, 



uniformly for t in compact sets of IR n . Hence, d Xj d Xk (p G D(e a( - X ^) for any a > 0. 

We next turn to third order derivatives. Consider the derivative of — Atp + (V — £)(p = 0. 
For any j G {1, n}, 

d Xj A(p = {d Xj V)<p + (V-S) d XjV . 

Since V is a polynomial, Proposition ()3.2|) shows that d x . Aip G _D(e a ^), for any a > 0. 
Thus, the function p i— > pj (X^=o Pf) ^(p) satisfies the Paley-Wiener conditions. 
Consider now any triple of indices j, k, I. For 52j=o Pj — R(t)-> we nave 



2 

r 

i=o 



| ^ + itj) (p k + it k ) (pi + iti) | < C(t) \pj + itj\ ^ P. 

3 

Hence, using this estimate with ()4.6|) . we deduce that 

/ I (Pj + %) {Pk + itk) ipi + Hi) | 2 \tp{p + it) | 2 dp 

jR™\B R{t) 

< C 2 (t) [ |Pi + %| 2 (X>?) \$<P + it)\ 2 dp 

J ]R n \B D(t\ \ „•_, / 



(*) \ i=l 



< OO, 



uniformly for t in compact sets of M n . Therefore, the Paley-Wiener Theorem asserts that 
d Xj d Xk d Xl (p G D(e a{x) ), for any a > 0. 

We now proceed by assuming D^ip G /}(e a ^), for any a > and any /?, such that 
|/3| < m. Let a have |a| = m + 1. Let <5 be any multi-index of length m — 1. Differentiating 
the eigenvalue equation again, Leibniz's formula yields 

D^Aip = C i (D^iV-S)) D^tp, (4.8) 

0<7<<5 

where the are multinomial coefficients. The induction hypothesis and the assumption 
that V is a polynomial show that D a Atp G L 2 (M n ). Therefore, p p a ($^j=i Pf) <p(p) 
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satisfies the Paley- Wiener conditions. In a, there are two indices, ctj and a^, not necessarily 
distinct, which are larger or equal to one, such that we can write 

(p + it) a (4.9) 

= (pi + iti) ai ■ ■ ■ (pj + itjp- 1 ■■■( Pk + it k ) ak - 1 ■■■(Pn + itn) an (p 3 + itj) (p k + it k ). 

Estimating the absolute value of the last two factors by ()4.7|) and using that 

a = (ai, ■ • • , atj — 1, • ■ ■ , a>k — 1, • • ■ , a n ) has length m — 1, we see that p a <p(p) satisfies 

the Paley-Wiener conditions. Hence, D a tp £ D(e a( - X ^) for any a > 0. I 
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